Spin dynamics in (III,Mn)V ferromagnetic semiconductors: the role of correlations 



00 
O 

o 

< 



I 

o 
o 



> 
O 

cn 

in 
o 
oo 
O 



M. D. Kapetanakis and I. E. Perakis 
Department of Physics, University of Crete, and Institute of Electronic Structure & Laser, 
Foundation for Research and Technology-Hellas, Heraklion, Crete, Greece 
(Dated: August 25, 2008) 

We address the role of correlations between spin and charge degrees of freedom on the dynamical 
properties of ferromagnetic systems governed by the magnetic exchange interaction between itiner- 
ant and localized spins. For this we introduce a general theory that treats quantum fluctuations 
beyond the Random Phase Approximation based on a correlation expansion of the Green's function 
equations of motion. We calculate the spin susceptibility, spin-wave excitation spectrum, and mag- 
netization precession damping. We find that correlations strongly afi'ect the magnitude and carrier 
concentration dependence of the spin stifi'ness and magnetization Gilbert damping. 
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Introduction — Semiconductors displaying carrier- 
induced ferromagnetic order, such as Mn-doped III-V 
semiconductors, manganites, chalcogenides, etc, have re- 
ceived a lot of attention due to their combined magnetic 
and semiconducting properties |^ . A strong response 
of their magnetic properties to carrier density tuning via 
light, electrical gates, or current [1,0,01 can lead to novel 
spintronics applications and multifunctional magnetic 
devices combining information processing and storage on 
a single chip. One of the challenges facing such magnetic 
devices concerns the speed of the basic processing unit, 
determined by the dynamics of the collective spin. 

Two key parameters characterize the spin dynam- 
ics in ferromagnets: the spin stiffness, D, and the 
Gilbert damping coefficient, a. D determines the long- 
wavelength spin-wave excitation energies, loq ~ DQ^, 
where Q is the momentum, and other magnetic prop- 
erties. D also sets an upper limit to the ferromagnetic 
transition temperature: Tj, cx I? [ij. So far, the of 
(GajMn)As has increased from K ^] to ~173 K 

[l|, |7|. It is important for potential room temperature 
ferromagnetism to consider the theoretical limits of Tc. 

The Gilbert coefficient, a, characterizes the damping of 
the magnetization precession described by the Landau- 
Lifshitz-Gilbert (LLG) equation [^. A microscopic 
expression can be obtained by relating the spin suscepti- 
bility of the LLG equation to the Green's function 
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average over a grand canonical ensemble and Sq = 
l/VNJ2j^j^~^^^^i where Sj are spins localized at N 
randomly distributed positions . The microscopic ori- 
gin of a is still not fully understood ^ . A mean-field cal- 
culation of the magnetization damping due to the inter- 
play between spin-spin interactions and carrier spin de- 
phasing was developed in Refs.(9l.[loj. The magnetization 
dynamics can be probed with, e.g., ferromagnetic res- 
onance [m and ultrafast magneto-optical pump-probe 
spectroscopy experiments [1, [13, 13,[l3|. The interpre- 



tation of such experiments requires a better theoretical 
understanding of dynamical magnetic properties. 

In this Letter we discuss the effects of spin-charge cor- 
relations, due to the p-d exchange coupling of local and 
itinerant spins, on the spin stiffness and Gilbert damp- 
ing coefficient. We describe quantum fluctuations beyond 
the Random Phase Approximation (RPA) [l^ [l^ with a 
correlation expansion [17| of higher Green's functions and 
a 1/S expansion of the spin self-energy. To 0(1/S'^), we 
obtain a strong enhancement, as compared to the RPA, 
of the spin stiffness and the magnetization damping and 
a different dependence on carrier concentration. 

Equations of motion — The magnetic properties can be 
described by the Hamiltonian [H H—HMp+Hcom where 
the mean field Hamiltonian HMF—^2i'k.n^^^^\in^^"' 

de- 
scribes valence holes created by aj^.^, where k is the mo- 
mentum, n is the band index, and ekn the band disper- 
sion in the presence of the mean field created by the mag- 
netic exchange interaction [l6| . The Mn impurities act as 
acceptors, creating a hole Fermi sea with concentration 
c/i, and also provide S = 5/2 local spins. 

H,,„ = (3cJ2AS^Ast^ + ^ ^(A^+Aslq + h.c.U2) 



where /? ~50-150meV nm'^ in (III,Mn)V semiconductors 
T| is the magnetic exchane interaction, c is the Mn spin 
concentration and Sq = l/V^Z^Wk "k+qn"kn' ^^'^ 
hole spin operator. AA — A ~ {A) describes the quan- 
tum fluctuations of A. The ground state and thermo- 
dynamic properties of (III,Mn)V semiconductors in the 
metallic regime [ch ~10^°cm~'^) are described to first 
approximation by the mean field virtual crystal approxi- 
mation, Hmf, justified for 5 ^ cxd [ij. Most sensitive to 
the quantum fluctuations induced by i?corr are the dy- 
namical properties. Refs. (ol. [isj treated quantum effects 
to 0{1/S) (RPA). Here we study correlations that first 
arise at 0(1/5^). By choosing the z-axis parallel to the 
ground state local spin S, we have = and = S. 



The mean hole spin, s, is antiparallel to S, s 
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The spin Green's function is given by the equation 

dt < S+ >= -2iS5{t) + /3c < (s X S_q)+ > 



/3c 



-^A « s+Q » +^ X 

< (cr„„' X ASp_k-Q)^ A[a[„ap„/] >, (3) 

kpnn' 

where A = j3cS is the mean field spin-flip energy gap 
and s = l/-/VX]kri '''nn/kn is the ground state hole spin, 
/kn = (ak„akn) is the hole population. The first line on 
the right hand side (rhs) describes the mean field pre- 
cession of the Mn spin around the mean hole spin. The 
second line on the rhs describes the RPA coupling to the 
itinerant hole spin flo|, while the last line is due to the 
correlations. The hole spin dynamics is described by 

{idt - ekri' + £k-Qn) < ajk-Qtaki > 
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When treating correlations in the realistic (III,Mn)V 
system, the numerical solution of the above closed sys- 
tem of equations of motion is complicated by the cou- 
pling of many momenta and bands and by unsettled is- 
sues regarding the role on the dynamical and magnetic 
anisotropy properties of impurity bands, strain, localized 
states, and sp-d hybridization [l|, |20|, |2l|, [22, 



23|. In the 

simpler RPA case, which neglects inelastic effects, a six- 
band effective mass approximation revealed an order 
of magnitude enhancement of D. The single-band RPA 
model jl5j also predicts maximum D at very small hole 
concentrations, while in the six-band model D increases 
and then saturates with hole doping. Here we illustrate 
the main qualitative features due to ubiquitous corre- 
lations important in different ferromagnets [3, US] by 
adopting the single-band Hamiltonian 15| . We then dis- 
cuss the role of the multi-band structure of (III,Mn)V 
semiconductors by using a heavy and light hole band 
model. 



In the case of two bands of spin-t and spin-| states 
[l5| . we obtain by Fourier transformation 
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The first term on the rhs gives the RPA contribution [13] , 
while the last two terms describe correlations. 

The correlation contributions to Eqs.© and ^ are 
determined by the dynamics of the interactions be- 
tween a carrier excitation and a local spin fluctuation. 
This dynamics is described by the Green's functions 
<C ASp_k-Q A[aj^^ap„'] whose equations of motion 
couple to higher Green's functions, ^ Sa^aa^a 3> and 
<C SSa^a describing dynamics of three elementary ex- 
citations. To truncate the inflnite hierarchy, we apply a 
correlation expansion and decompose <C Sa^aa)a ^ 
into all possible products of the form {a^aa'a) <^ S 
{S){aU) < flta >, (at a) < ASA[a^a] >, and (S) < 
a'aa^a where ^ a^'aa'^a is obtained after sub- 
tracting all uncorrelated contributions, (a^a) <C a^a 
from <§C a^aa^a ^ (we include all permutations of mo- 
mentum and band indices) • Similarly, we decompose 
<C SSa^'^a » into products of the form {SS) <^ a'' a 
{S){a^a) < 5 >, (S) < ASA[a^a] >, and (a^a) < 
ASAS" This corresponds to decomposing all opera- 
tors A into average and quantum fluctuation parts and 
neglecting products of three fluctuations. We thus de- 
scribe all correlations between any two spin and charge 
excitations and neglect correlations among three or more 
elementary excitations (which contribute to 0{1/S^)) 
In the case of ferromagnetic /3, as in the mangan- 
ites, we recover the variational results of Ref. [lH and 
thus obtain very good agreement with exact diagonaliza- 
tion results while reproducing exactly solvable limits (one 
electron, half filling, and atomic limits, see Refs.jisl.flit). 
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UJ + 6 + Srpa(Q, Uj) + Ecorr(Q, w) 



(5) 



where S — Pes gives the energy splitting of the local spin 
levels. Srpa is the RPA self energy 
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is the correlated contribution, where 



^ « S+A[ala,] » ^ ^ « AS^a\a^ » 



F '^lO-lOO meV is the hole spin dephasing rate I2g|. Sim- 
ilar to Ref.[l3| and the Lindblad method calculation of 
Ref.Q, we describe such elastic effects by substituting 
the spin-flip excitation energy A by A-l-iF. We obtained 
G and F by solving the corresponding equations to lowest 
order in 1/S, with pS kept constant, which gives E'^"" to 
0(1/ S'^). More details will be presented elsewhere [l8| . 

Results — First we study the spin stiffness D — D^^^- 

^corr _j_ ^corr 
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The RPA contribution D^^^ reproduces 
The correlated cotributions D'^" > and 
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FIG. 1: (Color online) Spin stiffness D as function of hole 
doping and interaction strength for the single-band model, 
c = lnm~^, r=0, Do = /2mhh, mnh = 0.5me. 



< were obtained to Oil/S"^) from Eq.® [ll: 
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where Q, k, and p denote the unit vectors. 

For ferromagnetic interaction, as in the manganites 
Til , 24 1, the Mn and carrier spins align in parallel. The 
Hartree-Fock is then the state of maximum spin and 
an exact eigenstate of the many-body Hamiltonian (Na- 
gaoka state). For anti-ferromagnetic /3, as in (III,Mn)V 
semiconductors, the ground state carrier spin is anti- 
parallel to the Mn spin and can increase via the scat- 
tering of a spin-| hole to an empty spin-| state (which 
decreases Sz by 1). Such quantum fluctuations give rise 
to D™" , Eq.jSl), which vanishes for /k^ = 0. 0"°" comes 
from magnon scattering accompanied by the creation of 
a Fermi sea pair. In the case of a spin-t Fermi sea, Eq.® 
recovers the results of Refs.[ll, Q. 

We evaluated Eqs.dl]) and Q for zero temperature 
after introducing an upper energy cutoff corresponding 
to the Debye momentum, = Btt^c, that ensures the 
correct number of magnetic ion degrees of freedom [l5| . 




FIG. 2: (Color online) Spin stiffness D for the parameters of 
Fig. [T] (a) and (b): two-band model, (c) and (d): dependence 
on the Fermi energy within the single-band model. 



Figs, mja) and (b) show the dependence of D on hole 
doping, characterized hy p = Ch/c, for two couplings (3, 
while Figs. [Tfc) and (d) show its dependence on f3 for 
two dopings p. Figure [1] also compares our full result, D, 
with D^^^ and D^^^ + 0"°". It is clear that the cor- 
relations beyond RPA have a pronounced effect on the 
spin stiffness, and therefore on Tc cx D [II, 0] and other 
magnetic properties. Similar to the manganites 3, 24 1, 
£)corr ^ g destabilizes the ferromagnetic phase. However, 
D'^" strongly enhances D as compared to D^^^ [l^ and 
also changes its p-dependence. 

The ferromagnetic order and T^. values observed in 
(III,Mn)V semiconductors cannot be explained with the 
single-band RPA approximation [l^, which predicts a 
small D that decreases with increasing p. Figure [1] 
shows that the correlations change these RPA results in 
a profound way. Even within the single-band model, 
the correlations strongly enhance D and change its p- 
dependence: D now increases with p. Within the RPA, 
such behavior can be obtained only by includingmultiple 
valence bands As discussed e.g. in Rcfs.jll, 0], the 

main bandstructure effects can be understood by con- 
sidering two bands of heavy (rnhh—0.5me) and light ( 
TO/;i=0.086me) holes. D is dominated and enhanced by 
the more dispersive light hole band. On the other hand, 
the heavily populated heavy hole states dominate the 
static properties and Ep. By adopting such a two-band 
model, we obtain the results of Figs. EJa) and (b). The 
main difference from Fig. [l]is the order of magnitude en- 
hancement of all contributions, due to mih/rnhh — 0.17. 



Importantly, the differences between D and remain 
strong. Regarding the upper limit of due to collective 
effects, we note from Ref.[7| that is is proportional to D 
and the mean field Mn spin. We thus expect an enhance- 
ment, as compared to the RPA result, comparable to the 
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Fig-© compares a with the RPA result as function of 
p. The correlations enhance a and lead to a nonlinear 
dependence on p, which suggests the possibility of con- 
trolling the magnetization relaxation by tuning the hole 
density. A nonlinear dependence of a on photoexcitation 
intensity was reported in Ref. [H (see also Refs.^HH). 

We conclude that spin-charge correlations play an im- 
portant role on the dynamical properties of ferromag- 
netic semiconductors. For quantitative statements, they 
must be addressed together with the bandstructure ef- 
fects particular to the individual systems. The correla- 
tions studied here should play an important role in the 
ultrafast magnetization dynamics observed with pump- 
probe magneto-optical spectroscopy 13, 13, 14, 21, 2a |. 
This work was supported by the EU STREP program 
HYSWITCH. 



FIG. 3: (Color online) Gilbert damping as function of hole 
doping for different interactions (3. c — Inm^'^, F — 20meV. 



difference between D and D^^^. 

The doping dependence of D mainly comes from its de- 
pendence on Ef, shown in Figs. Wl,c) and (d), which dif- 
fers strongly from the RPA result. Even though the two 
band model captures these differences, it fails to describe 
accurately the dependence of Ep on p, determined by the 
successive population of multiple anisotropic bands. Fur- 
thermore, the spin-orbit interaction reduces the hole spin 
matrix elements 
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For example, |cr^„/| is maximum 
when the band states are also spin eigenstates. The spin- 
orbit interaction mixes the spin-]" and spin- J, states and 
reduces |ct^„/P. From Eq.© we see that the deviations 
from the mean field result are determined by the coupling 
to the Green's functions <C cr+„, A[at a„'] > (RPA), < 
A5"^CT+„, A[at a„'] » (correction to RPA due to S^' fluctu- 
ations leading to 1?^}°" > 0), and < AS'+crf„j, A[a,\a„'] > 
(magnon-Fermi sea pair scattering leading to D'^" < 0). 
Both the RPA and the correlation contribution arising 
from AS^ are proportional to (T^„/. Our main result, i.e. 
the relative importance of the correlation as compared to 
the RPA contribution, should thus also hold in the real- 
istic system. The full solution will be pursued elsewhere. 

We now turn to the Gilbert damping coefficient, a = 
2S/lj X Im < at w ^ [9]. We obtain to 

0(1/5*^) that a = a^^^ + a™", where a^^^ recovers 
the mean-field result of Refs d, [13] and predicts a linear 
dependence on the hole doping p, while 
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arises from the carrier spin-flip quantum fluctuations. 
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